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1. Introduction 



The purpose of this article is to discuss a simple linear algebraic tool which has proved itself 
very useful in the mathematical study of spectral problems arising in elecromagnetism and quan- 
tum mechanics. Roughly speaking it amounts to replacing an operator of interest by a suitably 
chosen invertible system of operators. 

That approach has a very long tradition and appears constantly under different names and 
guises in many works of pure and applied mathematics. Our purpose here is not to provide a 
historical survey but to present an account of a specific approach from a personal perspective 
of the authors. On one hand we hope to provide a source of systematic references for the 
practitioners of our type of spectral theory and, hopefully, to convince others of the usefulness 
of this method. We do not know, but find very interesting, if the method which has proved itself 
so successful in theoretical studies has a chance of being useful numerically. 

The key elementary observation goes back - at least - to Schur and his complement formula: 
if for matrices 

P R- \ f E E+ ' 



R+ ) \ 
then P is invertible if and only if E i- is invertible and 

(1.1) P-^ = E - E+EzIE- , EZl = ~R-P~'^R+ . 

In fact the equivalence of invertibilities of P and E |_ holds for systems with a non zero lower 

right hand corner (see Lemma ITT)) but since here we always start with P and choose R± we can 
normally consider these simpler systems. Sometimes, in the context of index theory one considers 

operators P which are never invertible. In that case the index of P is equal the index of E |_ 

which is trivial to compute if E |_ is a matrix - see H2.4I 

In the study of linear partial differential equations the use of enlarged systems appeared in 
Grushin's work [7] on hypoelliptic operators. In a different context they were used in the thesis 
of the first author and the ± notation comes from there - see ^2.21 for an explanation in the 
context of linear algebra. As is seen there it is essential that the system is non-self-adjoint For 
that historical, if somewhat personal reason, we refer to the problem 

J Pu + R^u^ = V 

(1.2) \ R+u = v+ 

P : Hi — » H2 , R- ■ H- — > H2 , R+ ■ Hi — > , 



^^That distinguishes it from the KKT (for Karush-Kuhn- Tucker) systems popular in numerical studies - see 
for instance 5^ - which seem to be related to il2.1l beIow. 
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as a Grushin problem. If it is invertible, we call it well posed and we write its inverse as follows 

u \^( E E+ \( V 
u- \ E, ; I V, 



(1.3) 



In this case we will refer to E |_ as the effective Hamiltonian of P. That effective Hamiltonian 

normally has its own physical interpretation as will be seen in examples in IHl^ and 15.41 

To illustrate this by a straighforward example consider an operator P : L^(M") 
defined as a convolution, Pu = K -ku, with K e L°°(R"). We can take H± = L^(R") and put 
R-U-{x) = — (27r)^"u_ (— x), the negative of the inverse Fourier transform, and R+u{^) = u(^). 

One easily checks that the resulting Grushin problem is well posed and that E y is given by 

multiplication by K. This of course is the effective Hamiltonian for the convolution operator 
which is invertible on if and only if g L°°. 

The main difficulty in constructing useful Grushin problems is the choice of suitable operators 
R± and of the spaces on which they act. As will be illustated below that depends on the situation 
even though one can notice some underlying principles. 

The paper is organized as follows. In Sj^we present in detail several simple examples showing 
different ways of constructing Grushin problems. In 321 we review basic linear algebra techniques 
which are useful when studying Grushin problems arising in spectral theory. In we also 
show a typical parameter dependent estimate. Trace formulae which are central in the study of 
classical/quantum correspondence are the subject of 21 we give the basic idea in the context 
of Grushin problems and use it to prove the Poisson summation formula, in a way which lends 
itself to many generalizations. Finally, in 13 we describe - without proofs - four advanced 
examples: a remark on Lidskii-Lusternik-Vishik perturbation theory for matrices |14j . the 
Peierls substitution of solid state physics (from the work of Helffer and the first author [H]), 
the quantum monodromy approach to the Gutzwiller trace formula, and the asymptotics of 
scattering poles in electromagnetic scattering by convex bodies (from earlier work of the authors 
|20|.|21p. It would be very hard to survey all the examples in which the Grushin problem appears 
explicitely - not to mention, those in which it appears implicitely - and we again made some 
personal choices. 

Acknowledgments. We would like to thank Steve Zelditch for suggesting a talk on Grushin 
problems during the semi-classical semester at MSRI: this paper is a direct result of that. We 
are also grateful to Michael Overton for the references to Lidskii's perturbation theory which are 
the basis of H5.ll 

The work of the second author was supported in part by the National Science Foundation 
under the grant DMS-0200732. He is also grateful to Universite de Paris-Nord, for its generous 
hospitality in October 2003. 

2. Simple examples 

We give five examples. The first two are purely linear algebraic and the third and fifth are 
intended to show how well known objects in mathematical physics fit in the Grushin problem set 
up. The fourth example relates Grushin problems to analytic Fredholm theory which is one of 
the basic tools of spectral theory. 

2.1. The Moore-Penrose pseudoinverse. If P : C" C™ is a linear transformation, its 
Moore-Penrose pseudoinverse is the unique transformation P+ : C™ C" satisfying 

PoP+ oP = P. P+ oPoP+ = P+ , 
(2-1) ^ ^ ^ 

(Po p+y = pop+ , {p+ o py ^ p+ o p . 
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If P has full rank then 

(P*P)-ip* n<m 
P*[PP*)-^ n>m ■ 

In general P^ can be expressed by using the standard singular value decomposition P — U*YiV, 
and inverting the non- zeros entries in S. It is closely related to least square problems - see ^ 
Lecture 11]. 

Another way to describe the pseudoinverse is as 



P^ 



kor(P)J 



Tim(P) ■ 



ker(P)^ im(P) is 



where Try is the orthogonal projection on the subspace V, since Pfkor(P) 
bijective. 

The pseudoinverse is a special case of E in H1.3(l . with Hi = C"" , H2 = C™, and for a natural 
choice of R±, related to the least squares method. Before describing it, let us give a general 
statement relating the Grushin problem to H2.1|) : 

Proposition 2.1. In the notation of H1.2|l and H1.3|l we always have 

EPE = E , 

and the following equivalences 

PEP = P ^ E-P = 
(PP)* = PE ^ [R-E-Y = R-E- 
{EPy =EP ^ {E+R+y = E+R+ . 

In particular, when the conditions on the left hold, E — P+, in the sense that equations in 1)2. l|l 
are satisfied. 

We can now choose R± so that the conditions in Proposition 12 . II are satisfied. For that we 
simply put 

= P{Hi)^ , H+ = ker(P) 

R^ : = P{Hi)^ ^ H2, R+: ker(P) = H+ . 

This is generalized in ij3.t)l in order to take into account small eigenvalues of {P*P)i and 
(PP*)i 

2.2. Non-self-adjoint eigenvalue problems. In the previous example we had i?^ = P_. For 

non self-adjoint problems that constitutes an unnatural restriction as shown by the following 
elementary example. 

Let J be the n x n upper triangular Jordan matrix: 



J,. 



1 for k = j + 1 
otherwise. 



Let 





























V ) 




I 1 / 
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Then Je+ = 0, J*e_ = 0, ||e±|| — 1, and we can set up the foUowing wcU posed Grushin problem 
for J - A: 



J -X R- 



c" e c — > c" 



R-U- = U-C- , R+u = (u, cj 



( ' ) 




( \ 


A 














A 


V ) 




I 1 / 



R+ 

One easily checks that E |_(A) = A", and that 

(2.2) E+V+ ^ v+e+{X) , E.v ^ {v, e.{X)) , e+(A) 



If we add a small matrix perturbation, eQ to J the same problem remains well posed and, 
using a Neumann series argument for matrices (see H3.4l for a similar argument), 

oo 

Et+{X) = E^+{X) + Y,{-l)'^'E^WQ{E{X)Q)''-'E+{X) , 

k=l 

with uniform convergence for |A| < < 1 and e < eo, for some ep > 0. Using (|2.2|) we consequently 
see that 



(2.3) 



Et+{X) = A" - e(Qe+(A), e_(A)) + 0{e') 



Hence when n is large and |A| < 1 there will be no spectrum near A for a generic perturbation 
Q. This is illustrated in FigureQ] The most dramatic perturbation is obtained by taking Q with 
a large inner product {Qe+, e_). 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Figure 1 . Eigenvalues of a small random perturbation of a 200 x 200 Jordan 
block matrix (blue), and of the perturbation Q — ee_ (8) e+ (red) 
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This example is a linear algebraic model of the first author's thesis ^7] where the — h notation 
was introduced. It was motivated by the sign in Hormander's commutator condition - see |18| 
and also |23 for a lighthearted introduction. It is reflected here by the fact that 

(2.4) [J, J*]e± = ±e± . 

This example will be revisited in a more general context in ^3.51 

2.3. Feshbach method. The Feshbach method which has been useful in the study of quantum 
resonances fits in the framework of Grushin problems discussed in this paper. To review it we 
follow [HI and we refer to that paper for pointers to the vast literature on the subject. 

Suppose that a Hilbert space H can be written as a direct sum H ~ TC TC , and that the 
quantum Hamiltonian under investigation decomposes under this splitting as 

Assume now that for z € f2, an open set in C, the operator (zl^^ — H^^) is invertible. Following 
[3] we define the resonance function 

which, in classical terminology reviewed in S^is just the Schur complement of zl^^ — H^^ in 
z-H. 

It then follows, by block Gaussian elimination, that for z ^ a{H^ 

zea{H) 0^ct(Gv(z)) 
and moreover it can be verified directly that 



TVV\ 



I (C-i/)-irfC = tr / acGv(C)Gv(C)~'dC, 7.(i) =^ + ee",0<t <27r, 



(2.5) tr 

-'72 "'to 

that is, the multiplicities agree. 

To see how the Schur complement, and hence also the Feshbach method, fit in the Grushin 
scheme we consider the following larger operator 

0^) ■■ -DQTC ^n®n\ i?+ = (F^ 0-^) , i?_ = Ji;^ 

If zF^ — H^^ is invertible then this problem is well posed and Gaussian elimination shows that 

E^+{z) = -(zF^ - H^^) + iJ"(zF^ - i?^^)-iiJ^^ = -G^{z) . 

The multiplicity formula follows from general principles described in ^3.1l but of course it is easy 
enough to verify directly. 

2.4. Analytic Predholm theory. Here we recall the discussion of the appendix in For the 
basic facts from functional analysis we refer to [S] for an in-depth treatment and to ^| Sect. 19.1] 
for a comprehensive introduction. 

A bounded operator P : Hi ^ H2 between two Banach spaces, is called a Fredholm operator 
if the kernel of P, 



and the cokernel of P, 



kerP =''{?ie Hi : Pu = 0} , 



cokeTP'^= H2/{Pue H2 : u e Hi} , 
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are finite dimensional. It then automatically follows (see for instance IT51 Lemma 19.1.1] or the 
comment after the proof of Proposition I2.2|l that PHi is closed. For Fredholm operators the 
index is defined as 

ind P — dim kcr P — dim coker P . 

We have the following 

Proposition 2.2. Suppose that that for some choice of R± the Grushin problem \1.2\ is well 

posed. Then P : Hi H2 is a Fredholm operator if and only if E ^ : iJ_ is a Fredholm 

operator, and 

ind P = ind E ^ . 

Proof. As for all well posed Grushin problems we have that E- are surjective, and Ej^, i?_ 
are injective. 

The equation Pu = v is equivalent to 

(2.6) u^Ev + E+V+ , = E^v + E^+v+ , 
where 1;+ = R-^-u. This means that 

E^ : im P > im E |_ , 

is surjective, and that it induces a bijective map 

eI : iJs/imP — > H^/imE^+ , . 

In fact, if E-v e im_E ^ then we use (j2.6(l to see that v E imP. 

On the other hand, 

E+ : ker£;_+ — *kerP 

is a bijection. In fact, if m € kerP then u = and E = and the map is onto, which 

is all we need to check as E+ is always injective. 
We conclude that 

(2.7) dim coker P = dim coker E y , dim kcr P = dim ker E |_ . 

In particular the indices are equal. □ 

For Fredholm operators we can always take H± to be finite dimensional: let = dim ker P 
and n_ = dim coker P and choose 

R- : C"- — ^H2, R+ : Hi — > C"+ , 

of maximal rank and such that 

P_ (C"- ) n im P = {0} , ker(P+ fkcr p) = {0} . 

In that case E |_ : C"" — > C"+ and its index is, of course, n+ — n_. This argument also shows 

that the index does not change under continuous Fredholm deformations of P, and that PHi is 
closed: by Banach's open mapping theorem the operators E, in H1.3|l (constructed using linear 
algebra only) are continuous. 

The following standard result is proved particularly nicely using the Grushin problem frame- 
work: 

Proposition 2.3. Suppose that for z G C C, a connected open set, A{z) is a family of 
Fredholm operators depending holomorphically on z. If A{zo)^^ exists at a point zq G Q.. Then 
9 z I— > A{z)~^is a meromorphic family of operators. 
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Proof. Let zq Cz Q and let V{zo) be a small neighbourhood of zq. We can then form a Grushin 
problem for P — A{zo) as described before the statement of the proposition. The same give 
a well posed Grushin problem for P = A{z) for z € ^(^o), if ^(^o) is sufficiently small. Since 
the index A{z) is equal to zero we see that n+ = n_ = n and E^^{z) is an n x n matrix with 
holomorphic coefficients. The invertibility of E^^{z) is equivalent to the invertibility of A{z). 

This shows that there exists a locally finite covering of il, {^j}, such that for z e flj, A{z) 
is invertible precisely when fj{z) ^ 0, where fj is holomorphic in fij. Since is connected and 
since A{zi) is invertible for at least one zi S shows that all fj's are not identically zero. 

That means that det E ^(z) is non- vanishing holomorphic function in V{zo) and consequently 

E i-iz)'-^ is a meromorphic family of matrices. Applying (|f .1(1 we conclude that 

A{z)-' = E{z) - E+{z)E^+{z)-'E^{z) 

is a meromorphic family of operators in V{zq), and since zq was arbitrary, in Q. □ 

2.5. Boundary value problems. Let P be an elliptic second order operator on a compact 
manifold, X, with an oricntable smooth boundary, dX. For the simplest example we could take 
P — ~&l + V{x) on [a, 5], in which case all the objects below are easily described. 
We want to pose a Grushin problem for the Dirichlet realization of P: 

Pdu = f e L^X), u\ax=0. 

We then put 

n = H^{X)r\H^{X), H- = Hi{dX), n+ = H^{dX). 
Let T : H^{dX) H'^{X) be an extension operator, with the following properties: 

Tv\qx^v, duTv\dx^O, 

where dv is the outward normal differentiation at dX. The operator T can, for instance, be 
obtained by introducing normal geodesic coordinates {x, y) in a collar neighbourhood of dX, 
y G dX, and putting 

Tv{x, y) = x{x) exp(a;2Aax)w(y) , -^dx > , 

where x G C^{[0, 5)), x = 1, near 0. 
We then define 

i?- : 7Y- — > Ji, , -R+ : Ji, — * 

(2-8) dcf dof „ 

K-U^ — ~Piu^, K-\-U — o,jU\qx ■ 

If we denote by Pjv the Neumann realization of P, 

Pat?/-/ (X), 9,Mrax=0, 

we have 

Proposition 2.4. Wii/i R± defined by (|2.8|l i/ie Grushin problem for Po is well posed when P^^ 
exists. The effective Hamiltonian is given by the Neumann-to- Dirichlet map: 

E^+ = N, N : a^utaxi— > u\9x , Pu = , 
where the existence of N is guaranteed by the invertibility of Pn ■ 

Proof. We can write (|1.3f) explicitely using the Green operator. Gat Pn^j ^^d the Poisson 
operator: 

QNf^u, Pu = 0, duu\Qx= f ■ 
It can be easily constructed from Pj^^. 
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Using this notation we have 

Ev^Gnv + T{{Gnv)\ox) 
E^v = {Gnv)\ox 

= Qnv+ + T o Nv-f. 
E-+V+ = {QNV+)\dx 

A direct verification proves the surjectivity. To prove injectivity we see that injectivity of Pm 
gives 

P(u-Tii_) = 0, d^{u-Tu^) = =^ u-Tu^ = 0. 
Since u\qx= this shows that u_ = Tu^lox^ 0, and hence u = 0, as well. □ 

A more systematic approach and one related to another use of two-by-two systems P].|12l 
Sect. 20. 4] can be described as follows. Suppose that P : C°°(X) — > C°°{X) is an elliptic operator 
of order m, and that we have two sets of boundary differential operators, with transversal orders 
< m, 

B, : C°°(X)^C°°(aX), j = ,J, 
Ck : C'^{X)^C'^{dX), k^l,--- ,K. 

For instance we can consider P = A, Biu = d;^u\gx, Ciu = u\gx, J ^ K = 1. 
We want to study the boundary problem 

(2.9) Pu = f in X, CkU = hk \ndX, k = l,--- ,K, 
assuming that the boundary problem 

(2.10) Pu = f inX, Bju = gk indX, ■■■,.], 

is well posed. To avoid technical issues involving Sobolev spaces (see 12', Chapter 20]) we will 
remain in the C°° category. We then put: 

Hi = c°°(x), rL2^c°°{x)®c°°{dx)^ , H--c°°(ax)^, K+ = c°°(9x)■^ 

writing 



( vx \ 

Vox 

V vj^x ) 



and define 



(2.11) 



f Pu \ 

ClU 



„ dcf 
R-U- ~ 



( ^ 



dct 



BlU 
Bju 



V Cku J 
We have the following formal 

Proposition 2.5. Suppose that the boundary value problem l|2.10|l is well posed. Then the 
Grushin problem 

Qu + R_u_ — V , R^u = w+ , 
obtained using the operators H2.11|l is well posed and the effective Hamiltonian, 

E_+ : C°°{dXy C°°{dXf 
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/ .1 ^ 








H 


\ Cku I 



is a generalization of the Neumann-to-Dirichlet map: 

(2.12) E + : \ : |^| : \ , Pu ^ , B,u = v{ , j ^ 1, ■ ■ ■ , J . 



For boundary value problems one of the basic issues is showing that, on suitably chosen 
spaces, the operator u i-^ {Pu, Ciu, ■ ■ ■ Cku) has the Fredholm property. By Proposition 12 . 21 that 
is equivalent to showing the Fredholm property of the operator (|2.12() . The reduction to the 
boundary described in Proposition 12 . 51 will furnish us with another example in 

3. Basic techniques 

Here we present some general results about systems arising from considering Grushin problems 
and examples showing how they can be used. We recall that a Grushin problem for an operator 
P : Hi ^ H2 is a system 



Pu + R^u^ = V 
R+u = v+ 

where i?_ : — > iJ, i?+ : H — > iJ+ ,. In matrix form we can write 

^ dcf / P R- \ , Hi®H^ ^ H2®H^ 



(3.1) 



, R+ ^ 

We say that the Grushin problem is well posed if we have the inverse 



that is 

C- ) ^ (1 ^! J ( ; 

In this case we will refer to E |_ as the effective Hamiltonian of P. 

3.1. Two by two systems. Here we consider an invertible system 

A=( > : Hi®H2 Hi®H2.. 



(3.3) , 

^" ) : Hi(SH2 ^ Hi(SH2. 

\ -D21 -D22 / 

We first recall the formula involving an expression known as the Schur complement in linear 
algebra and as the Feshbach operator in mathematical physics: 

Lemma 3.1. Suppose that A22 is invertible. Then Bn is invertible, and 

(3.4) = An - Ai2A^^^A2i . 

Proof. Using BnAn + B12A21 = / and B11A12 = —-812^22 we see that 

BiiAii - 511^12^22^^21 = / - {B12 + Bll^l2^22^)^21 = / , 

and the left inverse property is derived similarly. □ 
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We now allow the entries of A to depend on a parameter, and denote differentiation with 
respect to that parameter hy A A. The next lemma explicitely shows that the traces of 
BiiB^^ and A22A22 differ by terms not involving any inverses. In the case of holomorphic 
dependence on the parameter that means that these traces differ by holomorphic terms which 
disappear in contour integration. Before stating this precisely let us recall some basic facts about 
trace class operators - see or |12l Sect. 19.1]. 

If Hj are infinite dimensional Hilbert spaces, the operator A : Hi ^ H2 '\s said to be of 
trace class if the self-adjoint operator [AA*)^ : H2 H2 has a discrete spectrum, and 
Sjli Mi < If ^ is of trace class, and Bi ^ Hi ^ Hi, B2 : H2 — + H2 are bounded operators 
then ABi and B2A are of trace class. 

li H = Hi = H2, A is of trace class we can define the trace of A as follows. Let {ej}J!Li be an 
orthonormal basis of H , then 

00 

tr A = ^~^(Aej, ej) h , 

and this definition is independent of the choice of a basis. 

Finally, if A : _ff ^ iJ is of trace class and B : H ^ H \s bounded then 

(3.5) tv[A, B] = Q, [A, B]= AB - BA. 

Lemma 3.2. Suppose that the operators Aij are of trace class. Then, when A22 is invertihle, 
we have 

(3.6) tr Bj~/Bn = tr ^22^22 - tr AB . 

Proof. This is a straighforward computation based on the formulae, AB = I, B = —BAB, cyclicity 
of the trace, and Lemma l3. II fwe note that all Bij, and in particular Bn, are of trace class). We 
obtain 

tr Bf/Sii = tr A22A22 +tiEi+ tr E2 , 
El = —AiiBii — A12B21 ■ Hi — > Hi , 
E2 — —A21B12 — A22B22 ■ H2 — > H2 , 

and 

trEi +trE2 = - tr AB = tr AB . 

□ 

The relevance of this discussion for Grushin problems (which in principle have i322 = 0) will 
become apparent in the next subsection. 

3.2. From one Grushin problem to another. Suppose that we have a well posed Grushin 
problem H1.2|l with the inverse given by (|1.3|) . 

We want to check if another Grushin problem is well posed: 



(3.7) 



Pu + R-U- = V 

R+u = -Cj 



The corresponding operator will be denoted hy V : Hi® iJ_ H2® H+ . If the inverse exists 
we will denote it by 5, with the corresponding notation for the entries. 
The simple answer is given in 
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Proposition 3.3. The Grushin problem (|3.7|) is well posed if and only if the following system 
of operators obtained from the solution (|1.3() of the well posed problem H1.2|l . 



-E-R- E-+ 













o) 




o) 






- R+E 


:) 


-E_ 










(3.8) g . 

is invertible, that is if and only if the matrix of operators has a two sided inverse. In that case 

(3.9) g-' = 

where 

Proof. In place of H3.7|l we can consider a larger system 

+ R-U_ = V 

= v+ , 
= v+ 

in which v, 5+, and u_ are given, and u, m_, and are unknown. We can solve (|3.7|l by putting 
M_ = 0. Using ()1.3|l with we can write 

u = Eiv - R-it) + E+V+ 

(3.10) . 

u_ — E_ {v — R-U-) + E , 

or, since R+u ~ 

R+E+V+ - R+ER^ii^ =v+- R+Ei 

E — E^R^u^ = 7i_ — E^v . 

which in turn can be rewritten as 

u- \ _ f v+ - R+Ev 
v+ J \ M_ — E^v 

Hence the invertibility of g implies that H3.7|) is well posed. In fact, we first obtain u_ by inverting 
g and then u by using the first equation in (|3.1U|) . When u = we see that 

from which the equivalence and (|3.9|l follow. □ 

We illustrate Proposition 13 . 31 with an example which is also the basis for ^4.21 below. Let us 
consider 

(3.11) P = P{z)=^ hD^~ z, x e §1 =^M/(27rZ). 

We formulate a Grushin problem as in Sect. 2] where it was motivated by JUI- For that we 
want to find R±{z) so that 

(3.12) Viz) 1^^' ( ^^-^l^ ^-^'^ ] : H\S') x C L^i^') x C, 
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is invertible. Rather than give the answer in a "deus ex machina" manner we follow our original 
reasoning. First, a boundary condition 

u(0) , 

is a natural choice. Then we can locally solve 

(P - z)u = 
R+u — V ' 

by putting 

u — I+{z)v — exp{izx/h)v , — e < x < 27r — 2e . 
This is the forward solution, and we can also define the backward one by 

u = — exp(izx/h)v , —2n + 2e < xe . 

The monodromy operator Af (z, ft,) : C — > C, can be defined by 
(3.13) I+{z)v{Tr) ^ I^{z)M{z, h)v{-K) , 

and we immediately see that 

M(z, h) — exp{2niz/h) . 

We use I±{z) and the point it to work with objects defined on §^ rather than on its cover: a 
more intuitive definition of M{z, h) can be given by looking at a value of the solution after going 
around the circle. 

Let X e C°°(§\ [0, 1]) have the properties 

X(x) = l, -e<x<7r + e, x(a;)=0, -tt + 2e < a; < -2e , 

and put 

E+{z) = xI+{z) + {l-x)I-{z). 

We see that 

{P-z)E+ = [P,x]I+{z)~[P,x]I-{z) - [P,x]-I+{z)~[P,x]-I-{z), 
where [P, x\- denotes the part of the commutator supported near tt. This can be simplified using 

{i/h){P - z)E+ + (i/h)[P, x]-I-{z){I - M{z, h)) = 0, 
which suggests putting 

R-{z)^{i/h)[P,x]^I-iz), 

so that the problem 



(P-z)m + P_(z)u_ =0 
P-l_(z)u = V 



has a solution: 



u = E^{z)v 
M_ = E \-{z)v ' 

with E |_(z) = / — M{z, h) . In fact, it is much more natural, and easier for full-blown microfocal 

generalizations, to consider a different R+{z) so that, with symmetry reminiscent of M2.2I 

(3.14) P_ (z)w_ = u_e_ (z) , R^{z)u ~ (u, e^{z)) , e±iz,x) = {i/h)[P,x]±{exp{i • z/h)){x) . 

One can show that with this choice of R±{z), (|3.12|) is invertible and then 

" E{z) E+{z) 
E^{z) E^+iz) 



V{z)-'^E{z)^ 
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where all the entries are holomorphic in z, and E^{z), E \-{z), are as above. The operator 

E is the effective Hamiltonian in the sense that its invcrtibility controls the existence of 

the resolvent: 

(3.15) (P - = E{z) - E+{z)E_+{z)-^E_{z) . 

The invcrtibility is independent of x with the properties described above. Hence we can move 
to a singular limit in the choice of x deform tt to 0. That means that we consider the 
following Grushin problem (with suitably modified spaces): 

{i/h){hDx ^ z)u{x) ~ 5q{x)u^ = v{x) 
u(0+) = v+ 

In fact, we can write 

u = E{z)v + E+{z)v+ , U- = E-{z)v + E- — V{z)v^ , 

where 



(3.16) 



We finally come to an application of Proposition 13.31 In (|3.14|) it would be nice to be able 
to take e_(z) = e+(z), that is to have a self-adjoint Grushin problem. That would also simplify 
matters in more complicated situations. Hence suppose that 

e±(z, x) = f{x) eyi^{ixz/h) . 

We then have to consider the invcrtibility of the matrix Q in Proposition 13. 31 - which here is an 
honest 2x2 matrix. A brief calculation shows that for z g E, C/ is equal to 

-A B 



E{z)v{x) = ]l[o.27r[ / exp(i(a; - y)z/h)v{y)dy , E+{z)v+ = v+ exp(ixz//i)]l[o,2,r[ 



2tv 



E^{z)v = — exp{2Triz/h) / exp{—iyz/h)v{y)dy, E ^(z) = 1 — exp(27riz//i) , 



/ f{x)f{y)dydx, B= / f{x)dx 



and we observe that \B\^ = A + A. 

Hence the condition for invcrtibility becomes 

Re(Ae-'^*^/'') ^0, 

and that will always be violated for some z G M. Hence we cannot have a well posed Grushin 
problem for all z S M with e_ = e+ in (|3.14() . 

3.3. Iterated problems. The Grushin problems can be iterated and this is particularly impor- 
tant when the intermediate Grushin problems are formal and only after one or more iterations 
we obtain a well posed problem. An example of a useful formal problem will be given in 
Before giving an example of that we start with the following simple 



Proposition 3.4. Suppose that (|1.2|) is well posed with the inverse given by H1.3|) . // 

H+®V- — >H^®V+ 



E^+ 
N+ 

is invertible, with the inverse 



F F+ 
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then the new Grushin problem 

P R_N_ \ 

N^R^ ) ^ -^H^V^, 
is well posed with the inverse given by 

I E — Ej^FE_ Ej^Fj^ 

Proof. We need to solve 

Pu + R-N^u^ = V 
N^R+u — v+ 

Putting N-U- — U-, and R+u — v+, we obtain 

Pu + R-U- = V 
R^u = w+ 

which is solved by taking 

u = Ev + 
{t_ — E^v + E 

Recalling the definitions of u_ and this becomes 

E^+ N- \ f v+ \ _ f -E^v 
N+ )[-u^ )^[ v+ 

Solving this using J- gives the lemma. □ 

We will mention one concrete example for which iterated Grushin problems are useful. In the 
notation of Proposition 12 . 51 consider for X an open set in R"+^, with a smooth boundary fl, and 
put Pu = Au, Biu = u\q, and Ciu = Vu\a, K ^ J = 1, where ^ is a vectorfield. If V is not 
everywhere transversal to Q we obtain the oblique derivative problem and the operator l|2.12|l is 
not elliptic and not self-adjoint. As in ^7] one can then construct a new Grushin problem for 
that operator using the structure of the set where V is not transversal to ft. A "baby" version 
of that type of problem was presented on the level of matrices in W2.2\ 



3.4. A Grushin approximation scheme. Let iJ be a Hilbert space, and Hq a finite dimen- 
sional subspace with an orthonormal basis {ej}jLi. Let us introduce 

R+ : H^C^ , i?_ = i?; : ^ H , 

given by 

N 

{R+u)j = (u, ej) , = ^ u_ j-ej . 

We want to consider the Grushin problem for the operator 

P = I-T, T : H ^ H. 

In many interesting situations we can reduce the study of a differential operator to the study of 
/ — T by factoring out an invertible term. 

The following lemma is related to the example presented in tl2.1l 
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Lemma 3.5. Let tt be the orthogonal projection on the span of ej 's. With the operators R± 
given above, and P = (1 — ttT) the problem (|1.2() is well posed, and the matrix (11.31) is given by 

(3.17) f EO El \ f 1-n i?_ 



E°_ E°_^ J \ R+{I + T{1- n)) R+TR^-1 
Proof. We observe that 

(3.18) R-R+ = vr , = Idciv , 7ri?_ = i?_ , R+tt ^ R+ , 

which leads to an immediate verification of (|3.17() . □ 

We can now consider the problem for 1 — T and we have 

Proposition 3.6. // ||(1 — 7r)T|| < 6 < 1 then the Grushin problem (|1.2|l with P ~ 1 — T and 
R± as in Lemm,a \H.!A is well posed, and the effective Hamiltonian has the following expansion: 

oo 

E.+ = El+ + ^ R+T{{1 - TT)TfR_ . 

k=l 

Proof. This is a typical Neumann series argument. Using Lemma 13.51 and writing I — T = 
I ~iiT - {I ~ tt)T we see that 

I -T R^ \ _ f I -ttT R^ \ f ( {I~n)T 

R+ )-[ R+ J y^'^"®''" ~ [ R+T{I - 7t)T 

where we used ()3.18|l to multiply 

1 - TT i?_ \ / (/ - 7r)T 

i?+(/ + r(l-7r)) R+TR^-l )\ 

Hence 

E E+ \^-^( {I-^)T \ V l-TT i?_ 
E^ E^+ ) ^ R+T{I -t:)T o) \ R+{I + T{1 - n)) R+TR^-l 

A _ / 
B Q ) ^ \ BA''-^ 

we immediately obtain the formula for E y. □ 

The difficulty with the approximation scheme described here is the need for an orthonormal 
basis. In practice that is rarely given in theoretical and, especially, numerical problems. To some 
extend that can be remedied as follows. 

Let T be an operator with the property replacing the smallness of (/ — Tr)T: we postulate that 
there exists a finite set, {ej}jLi, with the following property: 

M 

(3.19) \fueH3t,eC,reH Tu^Y.^,e,+r, ||r|| < , -||ru|l < || f < C||ru|| . 

As before we would like to construct a well posed (in the sense that its stability constant is 
controlled, not just that it is invertible) Grushin problem for I — T . 

First, we need to modify the spanning set, {cj^^L^. For that we introduce the Grammian 
matrix, 

G = ((e;, ej))i<ij-<M ■ 



Since 
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It is positive semi-definite and hence can be diagonahzed. We then can, after a unitary (in C*^) 
"reorganization", assume that {cj} satisfy 

(ci, gj) = SijXj , Ai > A2 > • • • > Am > . 

Suppose that A^- > (e/C)^ for j < L. The condition number, ||G||||G~^||, of the Grammian, 
G, for {ej}j^i is now bounded by G^max|Aj|/e^ so we can use G, and its inverse, to form a 
well posed Grushin problem. For that we change ej to e^j ^J\j, and denote by tt the orthogonal 
projection onto the span of e/s. We easily see the following 

Lemma 3.7. Condition H3.19|l implies 

\\il-n)T\\<5 + e\\T\\. 

Proof. In the notation of (|3.19|) we write 

N 

j=M+l 

and 

M / AT \ 2 

\\Tu~J2^,e,\\^M<S\\u\\+iY,X,\t,\^\ . 

j=l \M+1 / 

Since Xj < e^/G^, and \\i\\i2 < G||rw|| the estimate follows. □ 

We can now proceed as in Proposition 13 . HI 

3.5. A typical estimate. Specific application of the Grushin problem scheme - see for instance 
|S1 - involve estimates, often depending on a parameter. We would like to illustrate this in a 
situation loosely related to the approximation scheme described in H3.4I and more concretely to 
the example in H2.2I 

Let us assume that P = P{h) : Ti ^ Ti is & bounded operator. Suppose that there exist two 
orthogonal projections ■n± = TT±(h) : TL TL satisfying 

(3.20) ||P*7r_||, IIPtt+II -0(/i), |K_P(/ - 7r+)|| - o(/i) , ||7r+P*(/ - ^_)|| = o(/i) , 

(3.21) ||P(/-^+)u|| >;i||(/-7r+)u||, ||P*(/-^_)«|1 >;i||(/-7r_)u||. 

Here by a{h) — o{h) we mean that lim/i^o a{h)/h = 0. 
We then have 

Proposition 3.8. With P and tt± with the properties described above we define 

n± =^m7r± , : H-^n, R+ : H ^ 'H+ , 

R*_R^=ld'H_, R-R*_=n^, R+Rl = ld-H+, P;P+ = 7r+ . 

Then for h small enough the Grushin problem 

Pu + R^u^ = V 
R^u = V-\. 

is well posed and 

(3.22) h\\u\\ + \\u-\\ < G(||i;|| +/i||i;+||), 
where G is independent of h. 
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Proof. We start by rewriting our Grushin problem as 

J Pu + R^u^ = V 

(3.23) I =0 

~ dcf 7-,^ ^ dcf ^ j-,^ 

u ~ u — iij^v^ , V = V — PKj^v^ . 

We observe that 7r_i?_ = i?_ and i?+7r+ = _R+. Taking the inner product of the first equation 
in H3.23|l with Pu and using {I — tt+)u — u in (|3.20() . gives 

i {h^\\u\\ + \\Pu\\^) < \\Pu\\^ < ||^'w||||w|| -Re(P*7r_i?_u_,M) 

(3.24) <^\\Pu\\^ + \\ir + o{h)\\u_\\njn 

<^\\Pur + 2M^ + oih^)M' + oil)\\u_\\l,_, 

Hence also 

Wu^Wn^ = WR^u^W < \\Pu\\ + M < Cll^ll +o(l)|l^._|l«_ , 

and consequently 

\\u^\\h- + h\\u\\ < C\\v\\ . 

To obtain ifX^ we estimate ||u|| using = 0{h)\\v+\\H+- 

\\v\\ < \\v\\ + 0{h)\\v+\\H^ . 

Since by definition ||7i|| < + the estimate follows. 

This shows the injectivity and to sec the surjectivity we apply the same proof to the adjoint 
Grushin problem, observing that the assumptions are symmetric. □ 

The estimate (|3.22|l is natural and appears under different assumptions (see for instance [211 
Lemma 5.2] for another elementary abstract estimate). In the proof we could have considered 
h — 1 since we can scale h out of the hypotheses: 

f h-^ \ f P{h) i?_ \ / 1 \ / h-^P{h) i?_ \ 
I )\R+ I )\{) h R+ I ) ■ 

3.6. Application to pseudospectral estimates. To see the estimates of ij3.5l in use we relate 
them to the example in >i2.2l The general phenomenon observed there is the growth of the 
resolvent of a non-normal operator away from the spectrum, and the consequent instability of 
eigenvalues - see [TH j .p ^ .pS ] . 

Thus consider a general n x n matrix A. Let us then put P = P{X) = A — X, and 

TT- = 7r-(A) =^ ]lp(A)p(A)*</i2 , 71'+ = vr+(A) =^ lp(A)*p(A)<h2 , 

where for a selfadjoint matrix i?, TiB<r denotes the orthogonal projection on the span of eigen- 
vectors of B with eigenvalues less that or equal to r. 
We then see that the hypothesis (|3.2()|l are satisfied: 

\\Pn+u\\^ = (P*F7r+M, vr+u) < h^\\TT+u\\ , \\P*tt_u\\^ = (PP*7r_w, 7r_w) < h^\\7r_u\\ , 

\\P{I - 7r+)uf = {P*P{I - 7r+)u, {I - 7r+)u) > h^\\ {I - 7r+)u\\ , 

\\P*{I -'K-)uf = {PP*{I -7r-)u,{I -Tr-)u) > h^\\{I -Tr-)u\\ , 

7r_P(/-7r+) = 0, 7r+P*(/- 7r_) = 0. 

To see the last identities we note that 

P : ker(P*F - r) > ker(PP* - r) , P* : ker(PP* - r) > ker(P*P - r) , , 
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SO that 

lpP*</l2-P]lp»P>;j2 =0, Ip* p<h2Ptpp,-^f^2 =0. 

This shows that we can apply Proposition 13 . 81 and consequently that the Grushin problem con- 
structed there has the inverse: 



E{X,h) E+{\h) 
E_{\h) E^+{\,h) 



\ f 0{l/h) 0{l) \ . j,„(A,/i) ^ (p„ „ (pn(A,/i) 

0{l) 0{h) j-^tB^ ^luw^ 
where n(X,h) = tT^\A-X)*(A-x)<h^- Using we see in particular that 
(3.25) \\{A - A)-i|| ^ \\E^+{X,h)-'\\ + Oil/h) , 

since ||i?+ti+|| ~ \\v+\\, and ||i?lM_|| ~ Hu-H, where a ^ h means that h/C < c < Cb for a 
constant independent of h. 

In the example presented in H2.2I where A was equal to a Jordan block matrix, we can take 
any |A|" <C /i < |A| to obtain a Grushin problem with n(A, h) = 1 and E |_(A, h) = A". 

4. Trace formulae 
4.1. Basic idea. Suppose that P = P(z). Writing dzA(z) = A{z) we have 

£{z)^-£{z)'P{z)£{z), 

which gives 

(4.1) E^{z)PE+{z) = -E^+{z) - E^{z)R^{z)E.+{z) - E_+{z)R+{z)E+{z) . 
We recall that, formally, 

P{z)-' = E{z) - E+{z)E^+{z)-^E^{z) . 
Hence, assuming that we have no difficulty in taking traces, we obtain 

(4.2) tvP{z)P{zy^ = tv E^+{z)E^+{zy^ + tv E^{z)R^{z) + tv R+{z)E+{z) + tr PE , 
which is a special case of Lemma [3. 21 This gives 

Proposition 4.1. Suppose that P — P{z) is a family of Fredholm operators depending holo- 
morphically on z Cz where 17 C C is a connected open set. Suppose also that the operators 
R± — R±{z) are of finite rank, depend holomorphically on z Cz fl, the corresponding Grushin 

problem in well posed for z £ f2, and that E |_(zo)^^ is invertible at some zq £ £7. Let g be 

holomorphic in f2. Then for any curve 7 homologous to in Q., and on which P{z)~^ exists 
P[z)P{z)^^g(z)dz is of trace class and we have 



(4.3) tr / P{z)P{z)-^g{z)dz :^tT / E_+{z)E_+{z)-^ g{z)dz . 

Proof. Since EZ]. is a finite matrix for z G 7 we have that 

P{z)P{zy^g{z)dz =^ - I P{z)E+{z)E^+{z)-^E^{z)g{z)dz 

is an operator of trace class, and, arguing as we did before the statement of the proposition we 
obtain (gSll. □ 
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The condition that E ^(z) is a finite matrix is often too restrictive. To illustrate this in 

a simple example we use the results of i j2.5l Let and Ajv be the Dirichlet and Neumann 
Laplacians on a bounded domain X, with a smooth boundary dX. We now put Pm{z) = — A, — z, 
• — D.N. As described in ij2.5l we have a well posed problem for Pd{z) if P^'^{z) exists, and 

in that case E = N{z), the Neumann to Dirichlet operator. Similary we have a well posed 

problem for Pn{z) if P^^{z) exists, and in that case E ^(z) = iV(z)^^, the Dirichlet to Neumann 

operator. Hence if is a contour homologous to in the region where PtvC^)^^ exists we get 
— tr /^^ Po{z)^^dz = tr /^^ N{z)N{z)~^dz. Strictly speaking we cannot apply Proposition 14. II 
directly but as N{z) is a Fredholm operator we can locally use an iterated problem with R± of 
finite rank. Our contour can be made a sum of contours involving only these local problems. 

Similarly we have j^^ PN{z)^^dz ~ 0. We can consider an analogous contour and write 
any contour 7 as + 7w- This leads to the following formula: 

(4.4) tr J {{-An - z)-^ - {-Ad - z)-^) dz ^ tr J N {z)-^ -^^N {z)dz , 

where N{z) is the Neumann to Dirichlet map for —A — z. For a non-trivial application of a 
similar idea in the context of resonances for the clastic Neumann problem see the work of Vodev 
and the first author |19| . 

4.2. Classical Poisson formula. To present an application of Proposition 14.11 we use it to 
derive the classical Poisson summation formula: 

(4.5) = E /(27rm) , / G C,- (K) , f{i) J f{x)e-^^^dx . 

Our proof here might well be the most complicated derivation of H4.5|l but as will be indicated 
in H5.2l it lends itself to far reaching generalizations. 

We start by rewriting (|4.5() using the operator P = hD^ on R/(27rZ): 

(4.6) tr f{P/h) J2 [ fi'/'^) (e'-^^"^) ' ^ (e^"^/'^) dz , 

The left hand side there can be written using the usual functional calculus based on Cauchy's 
formula: 

(4.7) tr/(^^^ = ^tr^/(^) (P-z)-idz, F = r+ - r_ , T±^^±iR, 

where we take the positive orientation of M and i? > is an arbitrary constant. We make an 
assumption on the support of the Fourier transform on /: 

(4.8) supp / C {-2ttN, 2nN) . 

We can now use the Grushin problem (|3.12() and its inverse given by H3.16|l . Applying Propo- 
sition ini with P{z) — {i/h){P — z) and g{z) — f(z/h) we obtain 

*'-K^)=-i/p^©*^'^^-^^^'^-^^^^"'^- 
We now use the expression for E |_ from ^3. 21 to write 

*'^(?) = 2^/ /(^)tra.M(z, /,,)(/ -M(z,/i))-idz 

+ f f(T)^'^9,M{z,h)M{z,hy\l-M{z,h)-Y'dz, 
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M{z,h) — cxp(27riz//i). The assumption H4.8|l and the Paley- Wiener theorem give 

\Kz/h)\ < e2"^|i'""l/''(Rez//i)-°°. 

Writing 

{I - M{z, h))-^ ^ J2 ^^(^' M'' + M{z, h)^{I - M{z, , 

fe=0 

for r_|_, and 

N 

M{z, h)-\l - A4{z, h)-^)-^ = ^ M(z, h)-^ + M(z, h)~^'^{I - M{z, h))'^ , 

k=l 

for r_, we can ehminate the last terms by deforming the contours to imaginary infinities (i? — > oo 
in H4.7|l '). and this gives H4.6() . 

4.3. An abstract version. In addition to demanding finite rank of R±, Proposition 14. II is re- 
strictive in the sense that we need to assume that the family of operators depends holoniorphically 
on the parameter z. Following 515! Appendix A] we present a result without that assumption. 

Let 7i be a complex Hilbert space and let us denote by £(7i, Ti.) bounded operators on Ti.. We 
consider §^ 3 1 1—^ A{t) £ £(7i, H), a closed curve of operators, in the sense of A{t) is strongly 
differentiable with a continuous derivative §^ 9 < i-^ A{t) e £{H,H). We write dA = Adt, and 
for another such 1 1— > B{t), 



B{t)A{t)dt = j BdAeC{n,n). 



If the values of A{t) are taken in an open subset V of £(7i,7i), we will will say that A{t) is 
contractible in V", if 9 i A{t) £ V has a extension D 9 z i-> A{z) e V, B = {\z\ < 1}, 
dB = §1. 

With this terminology we have 

Proposition 4.2. Suppose that 

(4.9) §1 9 t ^ V{t) = ^^^j^^^ ^ : H ® 7^ ® , 

is contractible in the set of invertible operators on Ti (B , with B) 3 z i-^ 'P{z) continuous with 
values in operators of trace class. 
If P{t)-'^ exists for all t e then 

(4.10) tr j p-'^dP = tr y EZ\dE-+ , 

where we use the standard Grushin problem notation for the inverse of (|4.9|l . 
Proof. Schur's formula 1)1. 1|) and the fact that dP is of trace class give 

tr p-^dP tr EdP - tr E+EzXE_dP . 
Using H4.1|l and the cyclicity of the trace we see that 

tr p-^dP = trEzXdE^+ + uj , w =^ tr{dR+E+) + tr{E^dR-) + tr{EdP) . 
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r{z) 



We will obtain H4.10|l when we show that lo, which is defined on the circle, extends to a closed 
form in the unit disc. Since t ^ V{t) is assumed to be contractible we can use the same notation 
for its extension to the unit disc: 

■ P{z) R^{z) 
R-{z) R+-{z) 

Lemma \\\.2\ then shows that lo is. & restriction to the unit circle of a one form defined in the unit 
disc: 

uj = tv{dR+E+) + tv{E^dR^) + tr{EdP) + tr(di?+_£;_+) . 
To compute dtu we note that d£ = —£dV£ and consequently 

-dE = EdPE + E+dR+E + EdR^E^ + E+dR+^E^ , 

- dE+ = EdPE+ + E+dR+E+ + EdR^E^+ + E+dR+^E^+ , 

- dE^ = E^dPE + E^+dR+E + E^dR-E^ + E^+dR+^E^ , 

- dE^+ = E^dPE+ + E^+dR+E+ + E^dR^E^+ + E^+dR+-E^+ . 
Hence, using the natural notation for operator valued differential forms, we obtain 

duj = tvdR+ A EdPE+ + tvdR+ A E+dR+E+ 

+ tTdR+ A EdR^E^+ + ivdR+ A E+dR+^E^+ 

- ii E^dPE A dR^ -ti E^+dR+E A dR^ 

- ti E^dR^E^ A dR^ - tvE^+dR+_E_ A dR- 

- tr EdPE AdP - tr E+dR+E A dP 

- tr EdR-E^ A dP - tr E+dR+^E^ A dP 

+ trdi?+_ A E^dPE+ + trdi?+_ A E^+dR^E+ 

+ trdi?+_ A E^dR^E+ + dR+- A E^+dR+^E^+ . 

For a differential 1-form, /x, with trace class we clearly have tr /i A /i = 0. That shows that in the 
expression for duj the 2nd, 7th, 9th, and 16th terms vanish. Cyclicity of the trace also shows that 
the terms in pairs: (1st, 10th), (3rd, 6th), (4th, 14th), (12th,13th), (8th,15th) cancel each other. 
Finally, we also have tr/xi A fi2 = — A /ii, and it follows that the 5th and 11th terms cancel 
each other. Thus dtu = completing the proof. □ 

It is quite possible that Proposition 14.21 follows from some general topogical facts. It is not 
clear what are the weakest assumptions on P and dP to guarantee that / PdP is of trace class. 
For a discussion of one case of a weaker assumption see ^| Appendix A] . 

5. Advanced examples 

5.1. Around Lidskii's perturbation theory for matrices. In i j2.2l the equation for the eigen- 
values of the perturbation of one Jordan block is easily derived from H2.3|l : 

A"-eQ„i + eO(A)+0(e2) = 0, |A| < 1 . 

and the solutions are 

= ei/«|Q,^^|i/«e(2'^»^+^'-sO"i)/" + o(ei/") , l<£<n. 

Here we consider n fixed and are interested in the e — > asymptotics. 

In this section we will show how the Grushin problem approach applies to the study of pertur- 
bation of matrices with arbitrary Jordan structure. We restrict ourselves to an example suggested 
by Michael Overton which according to him contains the essential elements of the general problem 
studied in and [T?>j . 
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Let Ji be the 



upper triangular Jordan bloc matrix. Wc then consider 



(5.1) 

that is 



A = 







—> i 


J n 


Onn 


Onk 




J n 


Onk 




Ofe„ 


Jk 



k < 



where O^p denotes the £ x p zero matrix. 

The Grushin problem for A is a straightforward modification of the one for J„ in 

R- : — >c"ec"ec^ r+ -. c"ec"®c'^ — > c^ 

We then obtain the eflfective Hamiltonian, E |_(A) for A — X: 

A" 



A" 






A'= 
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and E±(X) are similarly constructed from the three e±(A) vectors. 
Suppose we now consider a perturbation of A: 

A,^A + eQ, Q = 

(5.2) 

Qtj . c". — > C"^ = n , i = 1, 2 , na = fc . 
As in ij2.2l we see that the effective Hamiltonian for the perturbation is 
EL+{X) = -B-+(A) - eE^{X)QE+{X) + 0{^) . 
The effective first order perturbation is easily checked to be 





1 


9 Til 


■^ril 


g.;.^i 


E_{X)QE+{X) = 




021 

-inl 


(022 








031 


gfi 





0{X), 



J 



pq 



where Q^j denotes the ij'th entry of the matrix Q 

Suppose that the matrix (Q^i)i<i j-,<2 is diagonalizable with eigenvalues qi and g2- Then the 
eigenvalues of are given by the values of A for which the following matrix is not invertible: 

/ A"-egi eQi3 \ 

\ eQsi eQ32 A'=-eQ33 / 
Since k < n, and both k and n are fixed, perturbation theory gives 

Proposition 5.1. The largest modulus eigenvalues of A^ for e small are given by 

A] = ei/"|gjf/"e(2"^+^'-g9.)/" + o(ei/"), l<£<n, j = l,2, 

where qj are the eigenvalues (assumed to be distinct) of the (Q^i)i<ij<2 part of the perturbation 
matrix in (|5.2|l . 

A finer perturbation theory for matrices of size given by the number of distinct Jordan blocks 
will (most likely) give the general results of [T^ and 
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5.2. Generalized Gutzwiller trace formula. Trace formula; provide one of the most elegant 
descriptions of the classical-quantum correspondence. One side of a formula is given by a trace of 
a quantum object, typically derived from a quantum Hamiltonian, and the other side is described 
in terms of closed orbits of the corresponding classical Hamiltonian. 

Here we follow |2()| and outline the structure of a formula which is derived using a formal 
Grushin problem. It is an intermediate trace formula in which the original trace is expressed in 
terms of traces of quantum monodromy operators directly related to the classical dynamics. The 
usual trace formulae follow and in addition this approach allows handling effective Hamiltonians, 
such as the one described in ^5.31 below. 

Let P be a semi-classical, self-andjoint, principal type operator, elliptic in the classical sense, 
with symbol p, and a compact characteristic variety, p~^(0). Let 7 C p~^{0) be a closed primitive 
orbit of the Hamilton flow of p. The simplest example, and one discussed in M.2\ P = hD^, on 
the circle, p~^{0) = {{x, 0) x E C T*§^, and the Hamilton vector field is dx- More interesting 
examples are P = —h^Ag — 1 on a compact Riemannian manifold, or P = —h^A + V{x) with a 
suitable F on M". 

We can define the monodromy operator, M{z, h) for P — z along 7, acting on functions in one 
dimension lower, that is, on functions on the transversal to 7 in the base. We then have 

Theorem 1. Suppose that there exists a neighbourhood ofj, 57, satisfying the condition 

(5.3) m e and exp tHp{m) = m , p{m) = , < |i| < TN =4> m g 7 , 

where T is the primitive period of'j. If f e C~(R), supp / C {-NT, NT) \{0}, x e C^{^), 
and A G ^'^^'"(X) is a microlocal cut-off to a sufficiently small neighbourhood ofj, then 

(5.4) tTfiP/h)xiP)A^— 5^ tr / f{z/h)M{z,h)''—Miz,h)xiz)dz + Oih^), 

— TV — 1 

where M(z, h) is the semi-classical monodromy operator associated to 7. 

The dynamical assumption on the operator means that in a neighbourhood of 7 there are no 
other closed orbits of period less than TN, on the energy surface p = 0. We avoid a neighbourhood 
of in the support of / to avoid the dependence on the microlocal cut-off A. 

The monodromy operator quantizes the Poincare map for 7 and its geometric analysis gives 
the now standard trace formulae of Selberg, Gutzwiller and Duistermaat-Guillemin. The term 
k = —I corresponds to the contributions from "not moving at all" and the other terms to 
contributions from going |fc -|- 1| times around 7, in the positive direction when fc > 0, and in 
the negative direction, when k < —I. For non-degenerate orbits the analysis of the traces on 
monodromy operators recovers the usual semi-classical trace formulae in our general setting - see 
pni Theorem 3]. 

The proof of the formula follows the lines of the proof of classical Poisson formula presented 
in M.2\ In the general situation where the circle is replaced by a closed trajectory of a real 
principal type operator we can proceed similarly but now microlocally in a neighbourhood of 
that closed orbit on an energy surfarce. The contour integral formula (|4.7|l is replaced by the 
Dynkin-Droste-Helffer-Sjostrand formula (see H, Chapter 8]) 

(5.5) tr / (^^ x{P)A = -\jj Q d.xiz){P ^T'A C{dz) , 

where x is a-n almost analytic extension of x, that is an extension satisfying dzxi^) — ^{\ Ini2;|°°) 
- see 1201 Sect. 6] and we want to proceed with a similar reduction to the effective Hamiltonian 
given in terms of a suitably defined monodromy operator. 
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To construct the monodromy operator we fix two different points on 7, toq, mi (corresponding 
to and tt in H3.12|l - (|3.13|l ). and their disjoint neighbourhoods, W+ and W- respectively. We 
then consider local kernels of P — z near mo and mi (that is, sets of distributions satisfying 
(P — z)u = near m^'s), ker„i^. (P — z), j ~ 0, 1, with elements microlocally defined in W±. and 
the forward and backward solutions: 

I±{z) : ker,„o(P-z) — > ker„^(P-z). 

We then define the quantum monodromy operator, M{z) by 

/_(z)M(z) = /+(z), Miz) : ker„,(P-z) ker,„„(P-z). 

The operator P is assumed to be self-adjoint with respect to some inner product (•,•), and we 
define the quantum flux norm on kermQ(P — z) as follows"'': let x be a microlocal cut-off function, 
with basic properties of the function x in the example. Roughly speaking x should supported 
near 7 and be equal to one near the part of 7 between W+ and W^. We denote by [P, x]w+ the 
part of the commutator supported in W-^ , and put 

{u,v)qf =^ {[{h/i)P,x]w+u,v) , u,v € ker„io(P- z) . 

As can be easily seen this norm is independent of the choice of x- This independence leads to 
the unitarity of Ai{z): 

{M{z)u,M{z)u) QF = (m, u)qf , u e ker,„(,(P - z) . 

For practical reasons we identify ker„ia (P—z) with I?'(R"~^), microlocally near (0, 0), and choose 
the idenfification so that the corresponding monodromy map is unitary (microlocally near (0,0) 
where (0,0) corresponds to the closed orbit intersecting a transversal identified with r*R"~^). 
This gives 

M(z,h) : I?'(R"-i) — > I?'(R"-i), 

microlocally defined near (0,0) and unitary there. This is the operator appearing in Theorem^ 
and it shares many properties with its simple version exp(27riz//i) appearing in (|3.13|) for 

As in i|3.2l we can construct a Grushin problem with the effective Hamiltonian given by 
E ^_(z,h) — I — M(z,h). However, now the problem is formal, that is all the inversion for- 
mulae are only valid microlocally* near 7. Since in Theorem^we are interested in taking traces, 
and not, for instance, locating eigenvalues or resonances, that is sufficient. 

Nevertheless, as one striking application of this point of view we can explain the way in which 
complex quasi- modes manifest themselves on compact manifolds a phenomenon which was 
already explicitely or implicitely noted in the works of Paul-Uribe, Guillemin, and Zelditch - see 
[221 and references given there. 

To explain it, let us recall the now classical fact (Lazutkin, Ralston, Colin de Verdiere, Popov) 
that for an elliptic closed geodesic on a compact manifold M one can construct approximate 
eigenfunctions concentrating on that trajectory, and that the corresponding approximate eigen- 
values are close to actual eigenvalues with arbitrary polynomial accuracy as energy increases. 
When the trajectory is hyperbolic that procedure no longer makes sense as the formal construc- 
tion of quasi-modes gives complex numbers. That can lead to the construction of resonances 
in scattering situations (Ikawa, Gerard, Sjostrand-Gerard) but cannot have a direct spectral in- 
terpretation when the manifold is compact. Despite that they make a direct appearance when 

iSee IIUI for an earlier mathematical development of this basic quantum mechanical idea. 
For a review of this important notion see |2ni Section 3]. Roughly speaking it corresponds to a localization 
of the behaviour of quantum states to relevant subsets of classical phase space. It does not guarantee global 
well-posedness in an honest Hilbert space sense. 
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traces are considered and we have the following consequence of recent work on inverse spectral 
problems (see |22| and [T^'l: 

Theorem 2. Let M be a compact Riemannian Riemannian manifold and 7 a closed hyperbolic 
trajectory of primitive length Lj . Let denote the sequence of eigenvalues of the Riemann- 
Beltrami operator, /ij. the sequence of complex quasi-modes associated to the trajectory 7, < 
Im^fc (well defined modulo 0(1 Re/ife|~°°)/ Suppose that for any m^'L \ {0}, mL-y is different 
from the length of any closed geodesic on M which is not an iterate of ^ . Then, for any m € Z\{0} 
there exists a neighbourhood Um of mL^ such that 

j k 

where both sums are meant in the sense of distributions on M, and e'''''* is defined only modulo 
C°°(R\0). 

In our approach, epecially in view of Grushin reductions to the effective Hamiltonians, it is 
important that we can consider operators with non-linear dependence on the spectral parameter. 
In that case, motivated by ProDOsition l4.1l the left hand side of H5.4|l is replaced by 

itr J f{z/h)d, [x{z) d,P{z) P(z)-i] AC{dz) , 

which for P{z) = P — z reduces to (|5.5|l . For a generalized version we refer to |2l)[ Theorem 2]. 

Finally we point out that the semi-classical Grushin problem point of view taken here, when 
translated to the special case of C°°-singularities/high energy regime, is close to that of Marvizi- 
Melrose and Popov (see references in [201) In those works the trace of the wave group was reduced 
to the study of a trace of an operator quantizing the Poincare map. 

5.3. Peierls substitution. In this section we will follow jQi to show how the Grushin reduction 
leads to a natural mathematical explanation of the celebrated Peierls substitution from solid 
state physics. It gives an effective Hamiltonian for a crystal in a magnetic field. For simplicity of 
the presentation we will consider the case of dimension two only, and of the first spectral band - 
we refer to [0] and JO] for the general case and for references to the vast literature on the subject. 

First we need to consider the case of no magnetic field. Mathematically this corresponds to 
considering a Schrodinger operator with a periodic potential: 

Po^-A + V, FeC°°(R2), V{x + j) = V{x), jer, 

where F is a lattice in M.^. In other words, 

(5.6) TaPo — PoTa , Tau{x) u{x — o) . 

The operator P is unitarily equivalent to a direct integral of Floquet operators, Pg, acting as P 
on Tig: 

He =^ {u G PL(K') : V7 G F u(a; - 7) = e''<«'^>w(x)} , 6 G R/F* , 
where F* is the dual lattice of F: 7* G F* 4=^ (7*, a) G 27rZ for all a G F. We denote by E and 
E* the fundamental domains of F and F* respectively. Explicitely, 



6PoC = j PgdO, 

{Bf){x, e)^Y. ^-'^'^'^f{^ - 7) , {C9){x) = -J— f g{x, 0)de , 
B : L^iR^) — > L^{R^/r*,He), C = B* = B'^ 
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The spectrum of Pq is absolutely continuous and equal to [Jf.^ff{^k{(^) ■ S e M^/F*}, where 
is the sequence of eigenvalues of Pg. Each interval in the union is referred to as a 
band and we assume that the first band is disjoint from all the other bands. 

We now want to find a Grushin problem for P — z which will be well posed near the first band. 
It turns out (see [HI Lemma 1.1]) that one can choose (j>{x,6), Pe(j){x,9) = Xi{6)(l){x,6), to be 
holomorphic, as a function of 9, in a complex neighbourhood of M"/r*. That implies that 

(koix) '^^^ iC(j)){x) , 

has very nice properties: |9"0o(a^)| < CqC"!^'/'-^. We then define the following Grushin problem: 
(5.7) 



Rl (i?«_)* , R°_u_{x) "-(7)0o(x - 7) . 

It is not hard to see that this problem is well posed for z close to the first band and away from 
all the other bands. The effective Hamiltonian is given by 

{E"_^{z)v+)ia) = J2{^S^,P'Eia-P))v+{f3), E{j) ^ -J— f \,{e)e:^'^'>) dO , 
per ^ ^ ' 

which is unitarily equivalent to the multiplication by z — Ai(0), the obvious effective Hamiltonian 
near the first band. 

The Grushin problem (|5.7I) does have the advantage of being stable under small perturbations 
and we will see it when the magnetic field is turned on. That correponds to adding a magnetic 
potential to our operator. Here we consider only a constant weak magnetic field B — hdxi Adx2- 

(5.8) Pb = {D,, ~ hx2f + Dl^ + V{x) , D,^ = ^9,^. , B = hdxi A dx2 ■ 

Although the operator Pb is no longer periodic in the sense of H5.6|l it commutes with magnetic 
translations: 

(5.9) T^Pb = PbTS, T^uix)":^' ei^^'-^^'>u{x-a), T^T^ ^ e-'^'''"''^^T^T^ . 
We now use the magnetic translations to modify the Grushin problem (|5.7|l : 

(5.10) (i?^^i_)(:r)''=^'^^._(7)rf0o(:^), (i??u)(7) ='(",rf<^o)L2(R.), 

Hl''^' {ueL\R^) : PBueL^R^)}. 
The operator Vb{z) commutes with 

^1 T-^„U^ l5^^P*(-B,aA7) 



r ^ 



7 



^,(a)'iii'e^^-"^"^^'u(a-7). 



It is shown in 9, Proposition 3.1] that when h is small {B = hdxi A dx2) then Vb{z) is invertible 
for z near the first band for Pq- Although it requires some technical work, roughly speaking it 
follows from the invertibility of and the smallness of the magnetic field. 
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The inverse has the same symmetries as Vb{z) and in particular t^E \-{B, z) = E z)t^ 

for all a G r. That implies that E B) is given by a "twisted convolution": 

(5.11) B)v+) (a) = e^^'''"^''^/B,.(« - P)v+{P) , |/b,.(7)I < Cel^l/^ . 

per 

Operators with kernels satisfying these properties form an algebra sometimes called the algebra 
of magnetic matrices. In 9j Proposition 5.1] it is shown that the inveribility of a magnetic matrix 
as an operator on £'^{T) is equivalent to its invertibility in the algebra of magnetic matrices. Let 
■Msif) denote the magnetic matrix associated to an exponentially decaying function on T, f: 

MB{f)ia, P) = et(s,"A/3)^(^ _ ^) ^ e*''^("^'3)^(^ _ ^) ^ 
where a is the standard symplcctic form on M^. It is easy to check that 

(5.12) MB{f)oMB{g)^MB{f #3 9), f#Bg{l)^ J2 e^<^'"^'^>/(a)5(/3)- 

Q+/3=7 

We are now getting close to the Peierls substitution which provides an elegant microlocal de- 
scription of E B). We can take the Fourier transform of an exponentially decaying function 

on r, /, 

/(0)1if^e'('^'^>/(7), 

to obtain a r*-periodic analytic function on R^. 

To simplify the presentation we assume now that F = Z'^. Then one can check 9, §6] the 
following fact: 

(5.13) OpTAf^^g) = OpTAf) ° Op^{g) , 

where Op^ denotes the semi-classical Weyl quantization of a function on M^: 

a(x,0 I — > a'"{x,hD) : L^{R) — > ^^(M) , 

provided that a and all of its derivatives are bounded (see 0). In view of H5.12|l and H5.13|l it is 
not surprising that the invertibility of Msif) in the algebra of magnetic matrices is equivalent 
to the invertibility of Op^(/) in the algebra of of pseudodifferential operators. This leads to 

Theorem 3. Suppose that the first spectral band of a Schrddinger operator with a 'L'^ -periodic 
smooth potential is separated from other bands, with 6 i— > E{6), the (2Tr7j)'^ -periodic first Floquet 
eigenvalue. Suppose that Pb is the corresponding magnetic Schodinger operator with B = hdxi A 
dx2. Then there exists {2'!rZ,)'^ -periodic (in 9) analytic function, E = E(9,z,h), such that for z 
in a neighbourhood of the first band, and h small 

zeaiPB) ^ Oea{OprAEi;z,h))), 

E{e, z, h) - E{e) -z + hEi{e) + h^E2{e, z) -f • • • . 

For the formulation for a general lattice and any dimension (in particular dimension three) we 
refer to |H1 where one can also find the discussion of the coefficients in the expansion of E{6, z, h). 
Considering the spectrum of the leading term, E{x,hDx), already shows how dramatic is the 
introduction of the magnetic field from the spectral point of view - see ^Hl and the references 
given there. 
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5.4. High frequency scattering by a convex obstacle. In this section we will outline the 
construction of a Grushin problem which reduces an exterior resonance problem to a problem 
on the surface of the obstacle. It was used in to describe the asymptotic distribution of 
resonances in scattering by a convex obstacle satisfying a natural (at least from the point of view 
of our Grushin problem) curvature pinching conditions. 

The study of resonances, or scattering poles, for convex bodies has a very long tradition going 
back to Watson's 1918 work on electromagnetic scattering by the earth. He was motivated by the 
description of the field in the deep shadow. It provided impetus for the work on the distribution 
of zeros of Hankel functions which are the resonances for the case of the sphere. For general 
convex obstacles the distribution of resonances was studied, among others, by Buslaev, Fock, 
Babich-Grigoreva, Bardos-Lebeau-Rauch, and Harge-Lebeau. We refer to for pointers to the 
literature on the subject. 

The problem can be described as follows. Let O C K." be a stricly convex compact set with a 
C°° boundary. We consider the Dirichlet (or Neumann) Laplacian on R" \ O, — ARnyo, and its 
resolvent, 

i?o(A) =^(-Ar„\o-A2)-i : L2(R"\C') — > i72(R» \ O) n iJi(R" \ O) , ImA>0. 

When we allow i?c)(A) to act on a smaller space with values in a larger space, it becomes 
meromorphic in A: 



Ro{X) : L2„„,p(R" \ O) ^ H^^iR" \ O) n Hl^^{R" \ O) , Ae 



C when n is odd 
A when n is even 



where A is the logarithmic plane. The poles of this meromorphic family of operators are called 
resonances or scattering poles. They constitute a natural replacement of discrete spectral data 
for problems on non-compact domains - see (24] for an introduction and references. 

The first step of the argument is a deformation of R" \ O to a totally real submanifold, F, with 
boundary dT = dO in C". The Laplacian — Ar,i\o on M" \ O can be considered as a restriction 
of the holomorphic Laplacian on C" and it in turn restricts to an operator on F, — Ar. When 
F is equal to e*^°M" near infinity then the resonances of — ARnyci coincide with the complex 
eigenvalues of Ar in a conic neighbourhood of R. That is the essence of the well known complex 
scaling method adapted to this setting. 

Normal geodesic coordinates are obtained by taking x' as coordinates on dO and x„ as the 
distance to dO . In these coordinates the Laplacian near the boundary is approximated by 

(5.14) Dl^ - 2xnQ{x', D^,) + R{x\ D^,) 

where R is the induced Laplacian on the boundary and the principal symbol of Q is the second 
fundamental form of the boundary. The complex deformation near the boundary can be obtained 
by rotating Xn in the complex plane: x„ ^-^ e^^Xn which changes (|5.14|l to 

(5.15) - 2e'''xnQix', D.,,) + R{x' , D,,) . 

The natural choice of 9 comes from the homogeneity of the equation: 9 = tt/3. 

It is also natural to work in the semi-classical setting, that is, to consider resonances of 
— /i^Arti^c) near a fixed point, say 1. Letting /i — > gives then asymptotic information about 
resonances of — ARn^Q. 

Hence we are lead to an operator which near the boundary is approximated by 

(5.16) Poih) = e-2W3((/j^^j2 ^ 2x„g(x', hD,,)) + R{x' , hD,,) , 
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and we are interested in its eigenvalues close to 1. Let us consider the principal symbol of (|5.1t)|) 
in the tangential variables. That gives 

Po{h) = e-2W3((/j^^j2 ^ 2xnQ{x', eO) + R{x',e) • 

We are interested in the invertibility of Po{h) — C for C close to 1 and that should be related to 
invertibility of the operator valued symbol pQ{h) — We rewrite it as 

Po{h)-C = hi (e-'^'/'{D'^+tt,) + X~z) , 

(5.17) ^ ^ 

t = h~ixn, X = h-i{R{x',e)~l), z^h-i{C~l), ^i^2Q{x',^'), 

that is, we rescale the variables using the natural homogeneity of Po{h) — (. On the symbolic 
level the operator (|5.1t)|) can be analyzed rather easily. We can describe {po{h) — C)~^ using the 
Airy function: 

(D^ +t)Ai{t) ^0, Ai{-Cj)^0, Ai e L^{[0,oo)). 

Thus we consider 

(5.18) Px-z = e-'^'''/^{D'^ + fit) + X-z, AgM, 1/C<h<C, |Imz|<Ci, 

where Ci will remain large but fixed. To simplify the notation we shall now put /i = 1 (all the 
estimates will clearly be uniform with respect to /z with all derivatives). 

Let > — Ci > — C2 > • ■ • > — Cfc > ■ • • be the zeros of the Airy function and let ej{t) = 
CjAi{t — Q) be the normalized eigenfunctions of 

e,(0) = 0. 

We recall that the eigenfuctions ej decay rapidly since for t — > +00 we have 

Ai{t) - (2V^)-4"3 exp(-2t5/3) . 
We now take N = N{Ci) the largest number such that 

|Ime-'^2V3^jv| <Ci. 
To set up the model Grushin problem we define 

< :i'([0,(^)) ^ C^, Rlu{j)^{u,e,), 1<J<N, 
^ ■ ' i?" :C^^L2([0,oo)), R'_^{Rlr. 

Using this we put 



^^■^^^ Bx = {ueLl : DIu , tue L^, u{0) = 0} , 

II"IIb.,a = Rcz)||w||i2 + \\DfU\\L2 + \\tu\\L2 . 

Since the eigenvalues of Px are given by A + e-^'^^/^Cj and are the corresponding eigenfunctions, 
we see that Vxi^) bijective with a bounded inverse. 

As in H5.3I our Grushin problem becomes "stable under perturbations". However, because of 
the rescaling, the symbol class of the inverse is very bad in the original coordinates: we lose 
when differentiating in the direction transversal to the hypersurface R—l — 0. Overcoming that 
requires some second microlocal techniques. Once that is in place the invertibility of Pa{h) — 
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(1 + hsz) for I Imz| < C is controlled by invertibility of an operator on the boundary with the 
principal symbol given by 

A-/i-^(i?(x',0-l), Ai-2Q(x',0- 

Here N depends on C which controls the range of Imz. 

The passage to a global operator on the boundary, E ^(2;), with poles of E \-{z)^^ corre- 
sponding to the rescaled resonances is rather delicate. We use [2., Section 6] a symbolic calculus 
which takes into account lower order terms near the boundary. This results in an effective Hamil- 

tonian, E i-{z), described in Theorem^] In a suitable sense it is close to the model operator 

described above. It has to be stressed that a restriction on the range of Re z has to be made: 

for every large constant L we construct a different E which works for |Rez| < L. The 

properties of the leading symbol remain unchanged but the lower order terms and the symbolic 
estimates depend on L. 

The detailed description of the effective Hamiltonian is quite technical and involves the second 
microfocal classes of pseudodifferential operators introduces in 1211 Section 4] . Nevertheless from 
a computational point of view the construction is quite straightforward relying on the Grushin 
problem described above and the Taylor expansion of the coefficients of the Laplacian (in normal 
geodesic coordinates) at the boundary. 

Theorem 4. Let W d (0,oo) be a fixed set. For every w S W and z £ C, \ Rez\ < 1/VS, 

|Imz| < Ci there exists E^, \-{z), a second microlocal pseudodifferential operator associated to 

= {p e T*dO : Rip) ^w}, N = N{Ci) such that for Q < h < ho{S): 

(i) If the multiplicity of the pole of the meromorphic continuation of (Agn^Q — C)"'^ is given by 
moiC) then 

(5.22) mo {h~'^ {w + hhj) = — tr (L E^ ^+{Sy^ —E,, ^+{S)dS , 0<e<l. 

27ri J\i-z\=e ' dz 

(a) If E^^ — <^T.^,h{Ew. vi'^)) [P'^h), p G T*dO, a-^^ji, the second microlocal symbol 

map, 

El^_^{z;p- h)^0{{\-Rez)). 
In addition for |A| < l/{C\/5) we have 

(5.23) ||i?°,_+(z;p;/^)-diag(z-A-e'2"/3C,(p))|U(c«,c«) <e«l, 
and 

(5.24) det£;° _+(z;p; h) = Q ^ z = X + e"^"/^0(p) for some l<j<N 
where the zero is simple. Here Qj{p) — Cj{'2Q{p))^ . 

(Hi) For |A| > l/{C^/5), ^ is invertible and 

El^_^{z-p- /i)-i=0((A-Rez)-i). 

In 1^ Section 9] we give a trace formula for E ^{z). For that we start with the obvious 

observation that the trace of the integral of E ^{z)~^{d/ dz)E |_(z) against a holomorphic func- 
tion / over a closed curve gives the sum of values of / at resonances enclosed by the curve. The 
proof of the trace formula involves a further Grushin reduction, a local lower modulus theorem 
and a good choice of contours. The gain is in obtaining an integral in the region where the 
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- Im J, =K^ 2 (ReX )3 -C 



Figure 2. The distribution of resonances for a convex obstacle satisfying the 
pinched curvature assumption (|5.25l) with jo = 1 . 



operator E is elhptic (roughly speaking in the pole free region). A good choice of /, yields 

an asymptotic formula (see |21l Theorem 1.2]) for the number of resonances resonances in bands 

kCj (Re A)^ - C < - Im A < KCjiRe X)^ + C , j < jo 

„_1 TT . 2 i ^1 

K = 2 3 cos — mm Q 3 , K = 2 3 cos — max Q 3 , 

where, as above, Q is the second fundamental form of dO and SdO the sphere bundle of 90, 
provided that we have the pinched curvature condition: 

(5.25) maxsaoQ^^Go+i 



TDlTLSdoQ V Oo 

Under this assumption the regions between the bands are resonance free - this is shown in Figure 
1^1 which illustrates the result. 
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